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ABSTRACT 

We have carried out a detailed study of one of the most popular evolutionary chan- 
nels for the production of Type la supernova (SN la) progenitors, the semi-degenerate 
channel (CO+MS), where a carbon/oxygen (CO) white dwarf (WD) accretes matter 
from an unevolved or slightly evolved non-degenerate star until it reaches the Chan- 
drasekhar mass limit. Employing Eggleton's stellar evolution code and adopting the 
prescription of Hachisu et al. (1999) for the accretion efficiency, we have carried out 
full binary evolution calculations for about 2300 close WD binary systems and mapped 
out the initial parameters in the orbital period - secondary mass (P-M?) plane (for 
a range of WD masses) which lead to a successful Type la supernova. We obtained 
accurate, analytical fitting formulae to describe this parameter range which can be 
used for binary population synthesis (BPS) studies. The contours in the P-Mi plane 
differ from those obtained by Hachisu et al. (1999) for low-mass CO WDs, which are 
more common than massive CO WDs. We show that white dwarfs with a mass as 
low as 0.67Af Q can accrete efficiently and reach the Chandrasekhar limit. We have 
implemented these results in a BPS study to obtain the birthrates for SNe la and the 
evolution of birthrates with time of SNe la for both a constant star formation rate 
and a single star burst. The birthrates are somewhat lower than (but comparable to) 
those inferred observationally. 
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1 INTRODUCTION 

Type la Supernovae (SNe la) appear to be good cosmo- 
logical distance indicators and have been applied success- 
fully in determining cosmological parameters (e.g. Q and 
A; Riess et al. 1998; Perlmutter et al. 1999). They are 
believed to be thermonuclear explosions of mass-accreting 
white dwarfs (WDs), although the exact nature of their 
progenitors has remained unclear. In the most widely ac- 
cepted model, the Chandrasekhar mass model, a carbon 
oxygen (CO) WD accretes mass until it reaches a mass 
~ 1.378M Q flNomoto, Thielemann fc Yokoi 1984| , close to 
the Chandrasekhar mass, and explodes as a SN la. One of 
the more popular scenarios in which this occurs is that a CO 
WD accretes mass by mass transfer from a binary compan- 
ion, which may be a main sequence (MS) star or a slightly 
evolved subgiant. This channel is referred to as the WD+MS 
channel (van den Heuvel et al. 1992; Rappaport, DiStefano 
& Smith 1994; Li & van den Heuvel 1997). Here the binary 
system originally consists of two MS stars, where the more 
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massive star evolves to become a CO WD by binary in- 
teractions. Subsequentely, when the companion has evolved 
sufficiently, it starts to fill its Roche lobe and to transfer 
hydrogen-rich materiai onto the WD. The hydrogen accreted 
is burned into helium, and then the helium is converted to 
carbon and oxygen. The CO WD increases its mass till the 
mass reaches ~ 1.378M0 when it explodes in a thermonu- 
clear supernova. Whether the WD can grow in mass depends 
crucially on the mass-transfer rate and the evolution of the 
mass-transfer rate with time. If it is too high, the system 
may enter into a common-envelope (CE) phase (Paczyhski 
1976), if it is too low, burning is unstable and leads to nova 
explosions where all the accreted matter is ejected. 

Hachisu et al. (1999, HKNU99), Hachisu, Kato & 
Nomoto (1999, HKN99) and Nomoto et al. fTMfl l have stud- 
ied the WD+MS channel for various metallicities. However, 
their approach was based on a simple analytical method 
for treating the binary interactions. It is well established 
(e.g. Langer et al. 2000) that such analytical prescriptions 
cannot describe certain mass-transfer phases, in particular 
those occurring on a thermal timescale, appropriately. Li & 
van den Heuvel 1199711 studied this channel (for a metallic- 
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ity Z = 0.02) with detailed binary evolution calculations, 
but only for two WD masses, 1.0 Mq and 1.2 Mq, while 
Langer et al. (j2HDII| investigated the channel (for metallici- 
ties Z — 0.02 and 0.001) considering Case A evolution only 
(where mass transfer occurs during core hydrogen burning) . 

The purpose of this paper is to study the binary chan- 
nel more comprehensively and to determine the detailed pa- 
rameter range in which this channel produces SNe la, which 
can be used for population synthesis studies. Employing the 
Eggleton stellar evolution code with the latest input physics, 
we construct a fine grid of binary models for a population I 
metallicity Z = 0.02 in sections 2 and 3, and then implement 
the results in a binary population synthesis (BPS) study in 
section 4. In section 5 we discuss the results and compare 
them to previous studies and in section 6 we summarize the 
main results. 



2 BINARY EVOLUTION CALCULATION 

In the WD+MS channel, the lobe-filling star is a MS star 
or a star not evolved far away from the MS. The star trans- 
fers some of its mass onto the WD, which grows in mass as 
a consequence. To determine whether the WD reaches the 
critical mass for a SN la, 1.378 Mq, one needs to perform 
detailed binary evolution calculations. Here we use Eggle- 
ton's (119711 119721 11973H stellar evolution code, modified to 
incorporate the WD accretion prescription of HKNU99 into 
the binary calculations. 

The Eggleton stellar evolution code has been updated 
with the latest input physics over the last 3 decades, as de- 
scribed by Han et al. itT^Hl and Pols et al. fE?931 IT59lfll . 
Roche lobe overflow (RLOF) is also treated within the 
code pan, Tout fc Eggleton 2000) . In our calculation, we 
use a typical Population I (Pop I) composition with hy- 
drogen abundance X — 0.70, helium abundance Y = 0.28 
and metallicity Z — 0.02. We set a = l/H p , the ratio 
of typical mixing length to the local pressure scaleheight, 
to 2, and set the convective overshooting parameter 5 ov 
to 0.12 {Schroder, Pols fc Eggleton T997] IPols et al. 19971 . 
which roughly corresponds to an overshooting length of 
~ 0.25 pressure scaleheights (-Hp). 

We adopt the prescription of HKNU99 for the growth 
of the mass of a CO WD by accretion of hydrogen-rich ma- 
terial from its companion. If the mass-transfer rate is above 
a critical rate, M cr , hydrogen burns steadily on the sur- 
face of the WD and, the hydrogen-rich matter is converted 
into helium at a rate M cr , while the unprocessed matter is 
lost from the system, presumably in the form of an opti- 
cally thick wind. The mass-loss rate of the wind is therefore 

Mwind = |M 2 | 

The critical mass-transfer rate is given by 



M c: 



5.3 x 10" 



(1-7 -X) 
X 



(Mwd-0.4), 



(1) 



where X is the hydrogen mass fraction and Mwd the 
mass of the WD (masses are in solar units and mass- 
accretion/transfer rates in M Q yr _1 ). If the mass-transfer 
rate is less than M cr but higher than M s t = \M cr , it is as- 
sumed that there is not mass loss and that hydrogen-shell 
burning is steady. If the mass-transfer rate is below iMcr 
but higher than ;M cr , hydrogen-shell burning is unstable, 




log P 1 (days) 

Figure 3. Regions in the initial orbital period - secondary mass 
plane (log P 1 , Mi) for WD binaries that produce SNe la for initial 
WD masses of 0.67, 0.70, 0.75, 0.80, 0.90, 1.0, 1.1 and 1.2 M Q . 
The region almost vanishes for = 0.67 Mq. 



triggering very weak shell flashes, where it is assumed that 
the processed mass can be retained. If the mass-transfer rate 
is lower than Mi ow = §Mcr, hydrogen shell flashes are so 
strong that no mass can be accumulated by the WD. There- 
fore the growth rate of the mass of the helium layer on top 
of the CO WD is 



M Ec = 77h|M 2 



where 




[Ma | > M cr , 

Mcr > |M 2 | > Mow, 

|M 2 | < Mow. 



(2) 



(3) 



When the mass of the helium layer reaches a certain value, 
helium ignites. If helium shell flashes occur, a part of the 
envelope mass is blown off. The mass accumulation efficiency 
for helium shell flashes according to HKNU99 is given by 



-0.175 (logM He + 5.35) 2 + 1.05, 
'/ii. = <; -7.3 < logA/Ho < -5.9, 

1 , -5.9 < logM He < -5. 



(4) 



The growth rate of the CO WD is then given by Mco, i.e. 



Mco = ?7h A4ho = nncVn\M2 



(5) 



We assume that the mass lost from the system carries away 
the same specific orbital angular momentum as the WD. 

These prescriptions have been incorporated into our 
stellar evolution code and allow us to follow the evolution of 
both the donor and the accreting CO WD. Altogether, we 
have calculated the evolution of 2298 CO WD binary sys- 
tems, thus obtaining a large, dense model grid. The initial 
masses, M\, of donor stars range from 1.5 to 4.0 Mq, the 
initial masses, M^u, of the CO WDs from 0.67 to 1.2 M , 
the initial orbital periods, P 1 , of the binaries from the min- 
imum period, at which a zero-age main-sequence (ZAMS) 
star would fill its Roche lobe, to ~ 30 d. 
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Figure 1. Examples of three representative binary evolution calculations. The solid, dashed, dash-dotted curves show the mass-transfer 
rate, M2, the mass-growth rate of the CO WD, Mco> the mass of the CO WD, Mwd, respectively, in panels (a), (c) and (e). The 
evolutionary tracks of the donor stars are shown as solid curves and the evolution of orbital period is shown as dash-dotted curves in 
panels (b), (d) and (f). Dotted vertical lines indicate the position where the WD is expected to explode in a SN la in all panels. The 
initial binary parameters and the parameters at the time of the SN la explosion are also given in each panel. 



3 BINARY EVOLUTION RESULTS 

In Figure^we present three representative examples of our 
binary evolution calculations. It shows the mass-transfer 
rate Mi, the growth rate of the CO WD, Mco, the mass 
of the CO WD, Mwd, the evolutionary track of the donor 
star in the Hertzsprung-Russel diagram (HRD) and the evo- 
lution of the orbital period. Panels (a) and (b) represent the 
evolution of a binary system with an initial mass of the 
donor star of = 2.00 M©, an initial mass of the CO 
WD of M WD = 0.75 M© and an initial orbital period of 
log(Pyday) = 0.20. The donor star fills its Roche lobe on 



the main sequence (MS) which results in Case A RLOF. 
The mass-transfer rate exceeds M cr soon after the onset of 
RLOF, leading to a wind phase, where part of the trans- 
ferred mass is blown off in an optically thick wind, while the 
rest is accumulated by the WD. When the mass-transfer 
rate drops below M a but it still higher than M st , the opti- 
cally thick wind stops and hydrogen shell burning is stable. 
The mass-transfer rate decreases further to below M at but 
remains above Mi ow , where hydrogen shell burning is unsta- 
ble, triggering very weak shell flashes, and the WD continues 
to grow in mass. When the mass reaches = 1.378 M©, 

the WD is assumed to explode as a SN la. At this point, 
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Figure 2. Final outcomes of the binary evolution in the initial orbital-period — secondary mass (log P 1 , AP,) plane of the CO+MS binary, 
where P 1 is the initial orbital period and AP, the initial mass of the donor star (for different initial WD masses as indicated in each 
panel). Filled squares indicate SN la explosions during an optically thick wind phase (IM2I > M CT ), filled circles SN la explosions after 
the wind phase, where hydrogen shell burning is stable (M cr > I-M2I > Mcr), filled triangles la explosions after the wind phase where 
hydrogen shell burning is mildly unstable (M st > IM2I > A^iow)- Open circles indicate systems that experience novae, preventing the 
WD from reaching 1.378 Mq, while crosses show systems that are unstable to dynamical mass transfer. 
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the mass of the donor is Af^ = 1.1356 Mq and the orbital 
period log(P SN /day) = -0.0975. 

Panels (c) and (d) of Figure show another example 
for an initial system with Mh, = 2.20 Mq, M^d = 0.80 M 
and log(Pyday) = 0.40. The binary evolves in a similar 
way as in the previous example and the binary parameters 
when the white dwarf reaches Mwd = 1.378 Mq are Mf N = 
0.8055 Mq and log(P SN /day) = 0.1614. The main difference 
between this example and the previous one is that RLOF 
begins in the Hertzsprung gap (i.e. after the secondary has 
left the main sequence; so-called early Case B evolution) and 
that, at the time of the explosion, the system is still in the 
stable hydrogen-burning phase after the optically thick-wind 
phase. 

Finally, the third example in panels (d) and (e) of 
Figure Q represents the case where mass transfer starts 
in the Hertzsprung gap and where the binary remains in 
the optically-thick phase even at the supernova stage. In 
this case, the initial binary parameters are M\ = 2.40 Mq, 
A4 d = 1.00 Mq and log(PVday) = 0.60. When M WD = 
1.378 M , they are M| N = 1.4834 M and log(P SN /day) = 
0.3610. 

Figure[3]shows the final outcome of the 2298 binary evo- 
lution calculations in the initial orbital period - secondary 
mass (logP 1 , M 2 ) plane. Filled symbols show that the evo- 
lution leads to a SN la, where the shape of the symbols 
indicate whether the WD explodes in the optically thick- 
wind phase (filled squares; (IM2I > M CI ), after the wind 
phase in the stable hydrogen-burning phase (filled circles: 
(M cr > IM2I > M CI ) or in the unstable hydrogen shell- 
burning phase (filled triangles; (M s t > > M\ ov ,). Sys- 
tems which experience nova explosions and never reach the 
Chandrasekhar limit and systems that may experience a CE 
phase are also indicated in the figure. 

In Figure and Figure |3 we also present the contours 
for the initial parameters for which a SN la results. The 
left boundaries of these contours in panels (3) to (8) are set 
by the condition that RLOF starts when the secondary is 
still unevolved (i.e. is on the ZAMS), while systems beyond 
the right boundary are systems that experience dynamically 
unstable mass transfer when mass transfer starts at the base 
of the red-giant branch (RGB). The upper boundaries are 
also mainly set by the condition of dynamical stability, since 
systems above the boundary have too large a mass ratio for 
stable mass transfer. The lower boundaries are caused by the 
constraints that the mass-transfer rate has to be high enough 
for the WD to be able to grow and that the donor mass is 
sufficiently massive so that enough mass can be transferred 
to the WD to reach the Chandrasekhar limit. SNe la are 
either produced through case A binary evolution (the left 
parts of the regions), in which hydrogen is burning in the 
centre of the donor star at the onset of RLOF, or early case 
B binary evolution (the right parts of the regions) , in which 
hydrogen is exhausted in the centre and the donor is in the 
Hertzsprung gap at the onset of the RLOF phase. 

In Figure |3 we overlay the contours for SN la produc- 
tion in the (log P 1 , M 2 ) plane for initial WD masses of 0.67, 
0.70, 0.75, 0.80, 0.90, 1.0, 1.1 and 1.2 M . Note that the en- 
closed region almost vanishes for Mwd = 0.67 Mq, which 
therefore sets the minimum WD mass for which this chan- 
nel can produce a supernova. In Appendix A, we present 
fitting formulae for the boundaries of these contours (accu- 



rate to 3 per cent for most boundaries), which can be used 
for population synthesis studies. If the initial parameters of 
a CO WD binary system is located in the SN la production 
region, a SN la is assumed to be the outcome of the binary 
evolution. 



4 BINARY POPULATION SYNTHESIS 

In order to estimate the supernova frequencies in the 
WD+MS channel, one first needs to determine the distri- 
bution and properties of CO WD binaries. There are two 
channels to produce such systems. In the first, the primary 
is in a relatively wide binary and fills its Roche lobe as an 
asymptotic giant, where mass transfer is dynamically unsta- 
ble. This leads to the formation of a common-envelope (CE) 
phase (Paczynski 1976) and the spiral-in of the core of the 
giant and the secondary inside this envelope (due to friction 
with the CE). If the orbital energy released in the orbital 
decay is able to eject the envelope, this produces a rather 
close binary consisting of the WD core of the primary (here, 
a CO WD) and the secondary. As is usually done in binary 
population synthesis (BPS) studies, we assume that the CE 
is ejected if 

QCEA£ orb > P b ind|, (6) 

where AP or b is the orbital energy released, Pbmd the bind- 
ing energy of the envelope, and «ce the common-envelope 
ejection efficiency, i.e. the fraction of the released orbital 
energy used to overcome the binding energy. We adopt 
-Ebind = Pg r — athPth, where Pgr is the gravitational binding 
energy, E t h is the thermal energy, and ath defines the frac- 
tion of the thermal energy contributing to the CE ejection. 

A second channel to produce close CO WD bi- 
nary systems involves so-called case BB binary evolution 
( Dclg ado fc Thomas 198l| l , where the primary fills its Roche 
lobe when it has a helium core, leaving a helium star at the 
end of the first RLOF phase. After it has exhausted the he- 
lium in the core, the helium star, now containing a CO core, 
expands and fills its Roche lobe again, transferring its re- 
maining helium-rich envelope, producing a CO WD binary 
in the process. 

The CO WD binary system continues to evolve, and 
the secondary will at some point also fill its Roche lobe; 
the WD will then start to accrete mass from the secondary 
and convert the accreted matter into CO. We assume that 
this ultimately produces a SN la if, at the beginning of this 
RLOF phase, the orbital period, Po r b, and secondary mass, 
M 2 , are in the appropriate regions in the (logP 1 , M2) plane 
(see Figure to produce a SN la. 

In order to investigate the birthrates of SNe la, we have 
performed a series of detailed Monte Carlo simulations with 
the latest version of the BPS code developed by Han et al. 
(2003). In each simulation, we follow the evolution of 100 
million sample binaries according to grids of stellar mod- 
els of metallicity Z — 0.02 and the evolution channels de- 
scribed above. We adopt the following input for the simula- 
tions (see Han et al. 1995 for details). (1) The star-formation 
rate (SFR) is taken to be constant over the last 15 Gyr or, 
alternatively, as a delta function, i.e. a single star burst. In 
the case of a constant SFR, we assume that a binary with its 
primary more massive than 0.8 Mq is formed annually. For 
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Figure 4. The evolution of mass-transfer rate, M2, mass-growth 
rate of the CO WD, McO; and the CO WD mass, JWwdj similar 
to panel (a) of Figure but using an analytical approach sim- 
ilar to the one of Hachisu et al. (1999) rather than full binary 
calculations. 




Figure 5. A comparison of the results of this paper with those 
obtained by Hachisu et al. (1999, HKNU99) and Hachisu, Kato 
and Nomoto (1999, HKN99). Dark contours show the parame- 
ter regions in the (logP, logMj) plane for different white-dwarf 
masses that lead to a SN la. The light grey contours are taken 
from HKN99 and HKNU99. 
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Figure 6. The evolution of birthrates of type la supcrnovae (SNe 
la) for a constant Pop I star-formation rate (3.5 Mg yr _1 ). Solid, 
dashed, dash-dotted curves are for a constant initial mass-ratio 
distribution and with «ce = a th — 1-0 (solid), 0.75 (dashed) and 
0.5 (dash-dotted), respectively. Dotted curves are with a^f; = 
a t h = 1.0 and for a mass-ratio distribution with uncorrelated 
component masses. Thick curves are based on the contours in 
this paper, while thin curves are based on the results of Hachisu 
et al. (1999). 
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Figure 7. Similar to Figure |5] but for a single star burst of 
1O 11 M . 



the case of a single star burst, we assume a burst produc- 
ing 10 11 M in stars. (2) The initial mass function (IMF) 
of Miller and Scalo RWIQl is adopted. (3) The mass-ratio 
distribution is taken to be constant or, alternatively, it is 
assumed that the component masses are uncorrelated. (4) 
We take the distribution of separations to be constant in 
log a for wide binaries, where a is the orbital separation. 
Our adopted distribution implies that ~ 50 % of stellar sys- 
tems are binary systems with orbital periods less than 100 

yr- 

The results of the simulations are plotted in Figures^ 

m 



5 DISCUSSION 

5.1 Comparisons with previous studies 

Langer et al. 1120008 considered the CO+MS channel con- 
sidering case A evolution for metallicities Z = 0.02 and 
Z = 0.001, respectively. They used a rather simple prescrip- 
tion for the WD mass accretion. We checked that we obtain 
very similar results, e.g. similar mass-transfer rates, similar 
final WD masses, if we adopt similar prescriptions. They also 
pointed out the possibility that binary systems with rather 
small initial white dwarf masses (~ 0.7 Mq) may produce 
SNe la. This agrees with our results that show that a binary 
system with A/w D = 0.67 Mq can produce an SN la. Li & 
van den Heuvel (19971 studied the evolution of WD binaries 
in search for progenitors of SNe la by using assumptions 
similar to those of Hachisu, Kato & Nomoto (1996; 1999) 
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Figure 8. The distribution of the initial masses of the CO WDs 
for the progenitors of SNe la according to the contours in this 
paper. The solid, dashed and dotted histograms are for SNe la 
produced at 0.3, 0.5 and 0.7 Gyr after a single star burst. The 
simulation uses a metallicity Z = 0.02, a constant initial mass- 
ratio distribution and Qce = Qth = 10. 



Roche lobe and the main-sequence star, respectively. In or- 
der to be able to make a comparison, we adopted a similar 
approach to calculate the evolution of a binary system with 
initial parameters M\ = 2.00 M , M^, D = 0.75 M and 
log(P'/day) = 0.20, the results of which are shown in Fig- 
ure 2] Comparison of Figure 0] and panel (a) of Figure Q 
shows that the estimate of the mass-transfer rate is signif- 
icantly different in this model compared to the results of 
our detailed binary calculations (also see Langer et al. 2000 
for a detailed analysis). The estimated mass-transfer rate 
is significantly larger, and the mass-transfer phase is cor- 
respondingly much shorter. The contours of HKNU99 are 
plotted as light grey curves in Figure Our contours are 
similar to theirs for massive WDs, but very different for low- 
mass WDs. For example, the contour for A/^d — 0.75 Mq in 
HKNU99 corresponds to significantly more massive donors 
(~ 2.7 M ) than our contour (~ 2.0 M Q ). The enclosed re- 
gion of the contour of HKNU99 is also larger. We consider 
our contours more physical as they are based on full binary 
evolution calculations rather than a simple prescription for 
estimating the mass-transfer rate and use the latest stellar 
evolution physics. 




0.7 0.8 0.9 1 1.1 1.2 



Figure 9. Similar to Figure but based on the contours of 
Hachisu et al. (1999). 

and Hachisu et al. (1999), but only for two WD masses, 
1.0 and 1.2 M . Our contours for the regions of WD bina- 
ries producing SNe la are consistent with theirs if we take 
into account that their model grid is much smaller and that 
their stellar evolution model parameters (e.g. overshooting 
parameter) may be different from ours. 

Hachisu et al. (1999, HKNU99), Hachisu, Kato and 
Nomoto (1999, HKN99) and Nomoto et al. itTMfll studied 
the WD+MS channel for progenitors of SNe la with differ- 
ent metallicities. They used the analytic fitting formulae of 
Tout et al. IU997H to calculate the radius and the luminosity 
of a main-sequence star. As the mass transfer proceeds on a 
thermal timescale for M^/M^d > 0.79, they approximated 
the mass-transfer rate as 

\M 2 \ = ^-M a x( CK \- CMS ,0) (7) 
tkh V Cms / 

where tkh is the Kelvin-Helmholtz timescale, and £rl and 
Cms are the mass-radius exponents of the inner critical 



5.2 Birthrates of SNe la 

Figure |S| shows Galactic birthrates of SNe la for the 
WD+MS channel. The simulations for a constant initial 
mass-ratio distribution give a Galactic birthrate of 0.6 — 
1.1 x 10~ 3 yr -1 , lower but not dramatically lower than 
those inferred observationally (i.e. within a factor of a few: 
3-4 x 10~ 3 yr -1 : van den Bergh & Tammann, 1991; Capel- 
laro & Turatto 1997). The simulation for the initial mass- 
ratio distribution with uncorrelated component masses gives 
a birthrate that is lower by an order of magnitude, as 
most of the donors in the WD+MS channel are not very 
massive which has the consequence that the white dwarfs 
cannot accrete enough mass to reach the Chandrasekhar 
limit. The birthrates based on our new contours (Fig.^J are 
0.6 — 0.8 x 10 -3 yr -1 , somewhat lower than the birthrates 
one would obtain whith the contours of HKNU99 for a con- 
stant mass-ratio distribution (0.8 — 1.1 x 10~ 3 yr -1 ). This 
follows directly from the fact that the parameter regions 
where low-mass WDs can produce a SN la is smaller in our 
study than in HKNU99 (see Fig.|KJ. 

Yungelson & Livio ( 199SJ gave a much lower SN la 
birthrate for the WD+MS channel, which is just 10% of the 
inferred SN la rate. HKN99 speculated that an important 
evolutionary path (case BB evolution) was not included in 
the model of Yungelson & Livio ljT998). However, Yungelson 
(private communication) has since pointed out that case BB 
evolution was indeed included though not mentioned espe- 
cially. Since we have also included this channel, we suspect 
that the differences between our result and theirs may be 
due to different treatments of CE evolution (see the dis- 
cussion of Han, Podsiadlowski, Eggleton 1995) and perhaps 
most significantly different assumptions concerning the ef- 
ficiency of hydrogen accumulation, where their efficiency is 
lower than ours. 

HKN99 proposed a wide symbiotic channel (i.e. 
WD+RG channel, in which a low-mass red giant is the mass 
donor) to increase the number of possible SN la progenitors 
and gave the corresponding regions in the (log P 1 , M2) plane. 
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However, at present it is not clear from a BPS point of view 
how to populate these parameter regions with WD systems. 
Our BPS model would predict that the Galactic SN la rate 
from the WD+RG channel is negligible (~2x 10 -5 yr -1 ). 
This is consistent with the conclusions of Yungelson & Livio 
(1998 ), although we think that this needs to be investigated 
further. 

Figure [7| displays the evolution of birthrates of SNe la 
for a single star burst of 10 11 Mq. Most of the supernova 
explosions occur between 0.1 and 1 Gyr after the burst. The 
figure also shows that a high common- envelope ejection effi- 
ciency leads to a systematically later explosion time, since a 
high CE ejection efficiency leads to wider WD binaries and 
it takes a longer time for the secondary to evolve to fill its 
Roche lobe. 

5.3 The distribution of the initial WD masses 

Figure 03 shows the distributions of the initial masses of the 
CO WDs that ultimately produce a SNe la according to our 
new calculations. The distributions are shown for systems 
that produce SNe la at 0.3, 0.5 and 0.7 Gyr after a single 
star burst. For clarity, we only show the distributions for 
the simulation with a constant initial mass-ratio distribu- 
tion and for «ce = «th = 1.0, as the other simulations gave 
similar results. The figure shows that there is clear trend 
for the peak in the initial WD mass distribution to move 
to lower masses with time, where the first SNe la preferen- 
tially have massive WD progenitors. Such a clear trend is 
not seen when the contours of HKNU99 are used (Fig. |j]J. 
Here the trend is not so clear, and the peak moves first to 
lower mass and subsequently to higher mass as the age in- 
creases. The difference between Figure Q3 and Figure H can 
be entirely understood from the different behaviour of the 
contours in the two studies and the fact that a massive donor 
in the WD+MS channel will evolve more quickly and hence 
produces a supernova at an earlier time. 

The brightness of a SN la is mainly determined by the 
mass of 56 Ni synthesized (Mm) during the explosion. As 
pointed out by Hoflich, Wheeler & Thielemann (1997) and 
Nomoto et al. (19915J, the amount of 56 Ni synthesized will 
generally depend on the C/O ratio in the core of the WD 
just before the explosion. Nomoto et al. (1999) speculated 
that the supernova brightness increases with the C/O ratio. 
We tested the C/O relation with mass in our stellar evolu- 
tion code and found, in agreement with earlier studies, that 
the C/O ratio decreases with increasing ZAMS mass and 
hence white dwarf mass. This would suggest that older pro- 
genitors, which tend to have lower initial WD mass (Fig. 01, 
should produce brighter SNe la in the WD+MS channel. 
However, irrespective of what the correlation of the C/O ra- 
tio with explosion energy may be, it is clear that the C/O 
ratio is important in determining the explosion energy. Since 
the C/O ratio depends on the age of the population and also 
the metallicity (see Hoflich et al. 1997, Langer et al. 2000 
for further discussions), this may introduce evolutionary ef- 
fects that may need to be taken into account in cosmological 
applications of SNe la. 

Nomoto et al. 1199911 have studied the distribution 
of the initial WD masses both for the WD+MS channel 
and the WD+RG channel. They found that the WD+RG 
channel tends to have higher WD masses on average as 



a more massive WD is generally needed to allow stable 
RLOF (rather than lead to dynamical mass transfer and 
a common-envelope phase). As a consequence, SNe la from 
the WD+RG channel would tend to be dimmer (assuming 
the above relation between supernova brightness and C/O 
ratio/WD mass). At an age of 10 Gyr, the WD+RG channel 
is the major single-degenerate channel to produce SNe la. 
This could explain why SNe la are much dimmer in elliptical 
galaxies. In our BPS simulations, the WD+RG channel was 
not important (see the earlier discussion) , and the WD+MS 
channel alone cannot explain the SNe la in elliptical galaxies. 
We will investigate this issue further, although at face value 
it might suggest that the theoretically less favoured double 



degenerate (DD) channel for SNe la (llben fc Tutukov 19841 
IWebbink fc Iben 1987) . in which two CO WDs with a to- 
tal mass larger than the Chandrasekhar mass limit coa- 



lesce, may be required (Yungelson & Livio 1998 


iHaii 1998 


|Yungelson & Livio 2000 


|Tutukov fc Yungelson 20021 . 


Finally, Langer et a 


1. (I2000H and Nomoto et al. (T1M?1 



pointed out that a more massive WD is needed to produce a 
SN la at low metallicity. As a low metallicity is in some sense 
similar to a young age from a stellar evolution point of view, 
this implies that the SNe la in low-metallicity populations 
should be dimmer. 



6 SUMMARY AND CONCLUSION 

Adopting the prescription of HKNU99 for the mass accre- 
tion of CO WDs, we have carried out detailed binary evolu- 
tion calculations for the progenitors of SNe la in the semi- 
degenerate channel (the WD+MS channel) and obtained the 
initial parameters in the (logP 1 , M2) plane that lead to SNe 
la. We provided fitting formulae that give the contours for 
the initial binary parameters that can be used in BPS stud- 
ies. We find that a CO WD with its mass as low as 0.67 Mq 
can accrete mass and reach the Chandrasekhar mass limit. 
By incorporating the contours into our BPS code, we ob- 
tained the birthrates of SNe la and their evolution with time. 
The Galactic birthrate is lower than but still comparable 
within a factor of a few to that inferred observationally. 
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APPENDIX A: FITTING FORMULAE FOR THE BOUNDARIES OF SN IA PRODUCTION REGIONS 

The upper boundaries of the contours in Figure [3] can be fitted as 

Mi = Mm(A + BlogP 1 + CtlogP 1 ) 2 + Z^logP') 3 , 3.2M WD ) (At) 
where 

A = -f .9 + 7.078Af{ VD - 2.284(M^ D ) 2 



B = -1.482 + 2.504Af WD - 0.5917(M^ D ^ 2 
C = 13.88 - 31.50M^ D + 16.86(M^ D 
D = 1.205 + 1.243M WD - 2.414(Mj 



(A2) 

\2 



In the fitting we ignored the vertical parts (at the right of the upper boundaries) as they are due to the finite spacing of our 
model grid. 

The lower and the right boundaries can be fitted as 

j _ ( 1.946 - 0.03261M^ D - 0.03639 log P 1 + 0.4366(log P 1 ) 2 if log P* < 0.24 

2 ~\ 1.145 - 0.007195M WD + 0.02563/(log P 1 ) 2 + 1.175 log P 1 + 0.7361(log P 1 ) 2 if log P i > 0.24 

The left boundaries are either defined by the onset of RLOF from the secondary on the ZAMS for high-mass WDs, or 
by the onset of RLOF at the beginning of the Hertzsprung gap for low-mass WDs. We also carried out additional binary 
evolution calculations and found that, when A/w D is lower than 0.73 Mq, SNe la cannot be produced by case A evolution. 
We therefore obtain 

. ; _ / -0.4096 + 0.05301M WD +0.02488M|, ifM WD > 0.73 

l0g ^ -\ 0.24, ifM WD <0.73 (A4) 

In the above fitting formulae, the masses are in Mq and the orbital period is in days. The formulae are valid for 
0.6 M < Mix D < 1.2 M , 1.5 M Q < M\ < 4.0 M and -0.5 < log(P i /day) < 1.5. If the initial parameters of a CO WD 
binary system is located in the SN la production region, i.e. the binary is below the upper boundary (equation IA1L above 
the lower and the right boundaries f equation IA3t and to the right of the left boundary ( equation I A4L a SN la is assumed to 
result from the binary evolution. The formulae are written into a FORTRAN code and the code can be supplied on request 
by contacting Z.H. 



